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I will shamelessly tell you what my bottom line is. It is placing 
balls into boxes. . . 

Gian-Carlo Rota, Indiscrete Thoughts 

Gian-Carlo Rota was a professor at MIT from 1959 until his death in 1999. 
He is arguably the father of the field today known as algebraic combinatorics. 
Rota had 50 students (notably Richard Stanley) and, as of this writing, he 
has 329 mathematical descendants, including many of the top names in the 
field today. (And me.) If putting balls in boxes was good enough for him, 
well, it should be good enough for any of us. 

In this article we will first give a gentle introduction to permutation statis- 
tics and Eulerian numbers, followed by an explanation of how to use the idea 
of balls in boxes in this setting. Then we will take those ideas and study a 
perfectly natural (though less well-known) refinement: the "two-sided" Eule- 
rian numbers. We will finish with a discussion of symmetries of the Eulerian 
numbers and a conjecture of Ira Gessel. 

Permutation statistics 

The symmetric group S n is the set of all permutations of length n, i.e., all 
bijections w : {1, 2, . . . , n} — > {1, 2, . . . , n}. We generally write a permutation 
in one-line notation: w = w(l)w(2) ■ ■ ■ w(n), so a typical element of Sj is 
w = 5624713. 

For any permutation w G S n , we define a descent to be a position r such 
that w(r) > w(r + 1), and we denote by des(w) the number of descents of 
w. For example, if w = 5624713, then there are descents in position 2 (since 
6 > 2) and in position 5 (since 7 > 1). Hence, des(w) = 2. The permutation 
12 • • -n is the only permutation with no descents, while its reversal, n ■ ■ ■ 21, 
has the maximal number of descents, with n — 1. 
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The function 

des: |JS n ->{0,l,2,3,...} 

n>l 

is an example of a permutation statistic. Any function from the set of permu- 
tations to a subset of the integers can fairly be called a permutation statistic. 
The reader is referred to [1] or [13] for more about permutation statistics gen- 
erally. 

Another common permutation statistic is inversion number, inv(w), which 
is the number of pairs r < s such that w(r) > w(s). Notice the useful fact 
that w and its inverse permutation must have the same number of inver- 
sions: if r < s and w(r) > w(s), then w(s) < w(r) and w~ 1 (w(s)) = s > 
r = w _1 (w(r)). For example, with w = 5624713, we get w' 1 = 6374125 
and we have inv(w) = inv(w _1 ) = 13. (Notice that des(u>) = 2 while 
des(io _1 ) = 3, so descent numbers do not enjoy this property.) Again, the 
permutation 12 ■ • • n has no inversions, while n • ■ • 21 has the most, with (™) 
of them. Here, every pair of indices is an inversion. 

We mention inv because we can use it to give an alternate definition 
of des as follows. First, recall the rth simple transposition, denoted ay, 
is the permutation that swaps r and r + 1 for some r and fixes all other 
elements of {1, 2, . . . , n}. There are n — 1 such elements in S n , one for each 
r = 1,2, . . . ,n — 1. The simple transpositions are very special because they 
form a minimal generating set for S n , i.e., every permutation w in S n can 
be obtained by applying a sequence of simple transpositions, and no proper 
subset of them will do the same. 

Now consider the permutation obtained by composing a simple transpo- 
sition with a permutation w: 

w o o r = w(l) ■ ■ ■ w(r — l)w(r + l)w(r)w(r + 2) ■ • • w(n). 

We see this action swaps the numbers in positions r and r + 1 of w . Thus r 
is a descent position of w if and only if r is not a descent position of wo ay. 
In particular, r is a descent of w if and only if inv(w) > inv(w o a r ), i.e., 

des(to) = \{r : inv(w o a r ) < inv(iu)}|. 

This identity is certainly not the most useful way to compute the descent 
number of a permutation, but it suggests that maybe we should consider the 
action of composition of w with a T on the left side as well. Here, we find 
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a r o w = (w 1 o a r ) 1 since a r 1 = a r , and so (remembering inverses have the 
same inversion number) we have: 

des(u? -1 ) = \{r : inv(w _1 o cr r ) < inv(itJ~ )}| 
= \{r : inv(a> o w) < inv(w)}|. 

Thus we may fairly call des(w) the number of right descents of w, while 
des(w~ 1 ) is the number of left descents of w. (The left descents are also 
called inverse descents in the literature, and this number is also denoted by 
ides(w). See, for example, |9] Section 13].) In this article we will study 
both the right descent number on its own, and how left and right descents 
are related. This leads, respectively, to the Eulerian numbers, and to the 
"two-sided" Eulerian numbers of the title. 



The Eulerian numbers 

With a little patience, or a bit of computer programming, one can compute 
the distribution of des over S n by calculating des(u>) for each w in S n . The 
reader may like to check by hand that in S3 there is one permutation with 
no descents, four permutations with one descent, and one permutation with 
two descents. With a little more effort, you can find the distribution for S4 
as well. As this task requires n\ things to check for a given n, you may not 
want to try anything bigger than n = 4 by hand. 

We write A n ^ for the number of permutations w in S n for which des(w) = 
i — 1. (This shift in indices is done for historical reasons; primarily it makes 
certain formulas come out nicer.) The numbers A n i are called the Eulerian 
numbers, after Leonhard Euler. In Tabled! we see the A n ^ for n < 8. 

There is a great deal of mathematical literature devoted to the distribu- 
tion of des, and we remark that this distribution manifests itself in many 
other ways. For instance, let runs(u>) denote the number of increasing runs 
of w (maximal increasing subwords), let asc(w) denote the number of ascents 
of w (positions r with w(r) < w(r + 1)), and let exc(u>) denote the number 
of excedances of w (positions r with w(r) > r). Then, remarkably, we have: 

A n<i = \{w e S n : des(w) = i - 
= \{w <E S n : runs(w) = i}\, 
= \{w E S n : slsc(w) — i — 1}|, 
= \{w G S n : ex.c(w) — i — 
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Table 1: The Eulerian numbers A. 



The second and third equalities are perhaps straightforward (there is nec- 
essarily one descent between each increasing run, and swapping each w(r) 
for n + 1 — w(r) trades descents for ascents), but that exc has the same 
distribution as des is not so obvious. 

Euler himself was not interested in permutation statistics, but was rather 
investigating solutions of certain functional equations in which the numbers 
An ti emerged. See the original [8j Caput VII, pp. 389-390] (available digitally 
online) or Carlitz's survey article, [I]. Carlitz points to several interesting 
occurrences of the Eulerian numbers in analysis and number theory in the 
first half of the 20th century, while in the early 1980s, Bjorner described the 
Eulerian numbers in a topological setting; see [2]. Suffice it to say, there are 
a variety of good reasons to study the Eulerian numbers, but for this article 
we will stick to the combinatorial point of view. 

We will now discuss one of the standard tools of enumerative combina- 
torics, with which we will encode the Eulerian numbers. Given a sequence 
of numbers, ao, a>i, a>2, ■ ■ we define the generating function for the sequence 
to be the power series (a polynomial, if the sequence is finite) in which the 
coefficient of f is 

a{t l = a + ait + a 2 t 2 + • • • . 

i>0 

So, for example, the geometric series 

= 1 + t + t 2 + ■ ■ ■ , 
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is the generating function for the sequence 1,1,1,.... Another example that 
might be familiar comes from the binomial theorem: 




so that for fixed n, (1 + t) n is the generating function for the binomial coeffi- 
cients ( n ) . Whether a generating function has finitely many terms or not, the 
reader should try not to worry about problems of convergence of the series. 
(For those who want to worry about it, rest assured that all series mentioned 
in this paper converge for \t\ < 1.) The reader should think of a generating 
function as a "clothesline" on which the sequence hangs. 

The generating function for Eulerian numbers, with n fixed, will be de- 
noted 

n 

A n (t) — An^t 1 = A n it + A n ^t 2 + • • • + A r ^ n t n . 
1=1 

We call this the nth Eulerian polynomial. For example, the fourth Eulerian 
polynomial is: 

A 4 (t) = t + lit 2 + lH 3 + t 4 . 

Notice that since A n j is the number of permutations with i descents, we can 
get the same generating function by summing not on i, but on the set of all 
permutations in S n , i.e., 

A n (t) = t des(w)+1 , 

w£S„ 

since the number of summands contributing t l to the sum is precisely the 
number of permutations with % — 1 descents. 

We will now discuss the "balls in boxes" approach to the Eulerian num- 
bers. 

Balls in boxes 

We begin with the generating function for the sequence a , ai, a 2 , . . ., where 
a k is the number of ways of putting n labeled balls into k distinct boxes. 

We claim the generating function for all the ways of putting n distinct 
balls into boxes is 

a + + a 2 t 2 H = ^ ^ '■ 

k>0 
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Indeed, if there are k different boxes, there are a total k n ways to put n 
labeled balls into the boxes; we have k choices for each ball. A typical choice 
with k = 9, n = 6 might look like this: 




We will use a shorthand notation for pictures like the one above, e.g., 

||56|2||14|||3 

Here we draw a vertical bar for the divisions between the boxes, so that the 
number of boxes is one more than the number of bars. Also notice that by 
way of standardization, we list the balls in each box in increasing order. So 
balls 5 and 6 were placed in the third box, ball 2 in the fourth box, and so 
on. We call such diagram for balls in boxes a barred permutation. 

We can partition the set of all barred permutations (arrangements of 
balls in boxes) according to the underlying permutations. Since we require 
numbers to increase within a box, we know that there must be a bar in each 
descent position of a barred permutation (i.e., at the end of every maximal 
increasing run), but otherwise, we can insert bars into the gaps between the 
numbers at will. So, for example, the barred permutations corresponding to 
w = 562143 are obtained from the barred permutation 

56|2|14|3 

by adding bars arbitrarily in the gaps between the numbers, i.e., there can 
be any number of bars to the left of 5, any number of bars between the 5 and 
the 6, at least one bar between the 6 and the 2, at least one bar between the 
2 and the 1, and so on. 

With the identity permutation, w — 12 • • • n, there are no required bars 
and we are free to add any number of bars in each of the n + 1 gaps be- 
tween the numbers. Thus, the generating function for barred permutations 
corresponding to 12 ■■■n is the generating function for the "multi-choose" 

numbers (( n ^ 1 )) = { k ^ n ) = (^t")' *- e -> ^ ne numDer of ways to choose, with 
repetition, k things from an (n + l)-element set. With the convention that 
each barred permutation has weight £ #ofbars+1 (= f# ofboxes ) ) we express this 
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generating function as: 

£ (" + r V = «- ,(i + t + f + -) „ d + t + e± 



^ — ^ bars in first gap bars in second gap 

= t ■ (1 + t + t 2 + ■ 



n + 1 gaps 

t 



(l_t)n+l 

For permutations with descents, we need only multiply this generating 
function by a power of t to reflect the number of bars required by descent 
positions: 



fc>0 v 



fc>o V n / 

^ + n - 1 - des(tw)V, ^ 



So, returning to the example of w = 562143, we get that the generating 
function for the sequence counting its corresponding barred permutations 
according to the number of bars is: 



t l . 



Summing (pQ) over all permutations in S n gives us, on the one hand, 

A n (t) 



(1 -t) n + 1 ' 

On the other hand, the sum gives the generating function for all barred 
permutations, which we know to be J2k>o k n t k - 
In other words, we have the following. 

Theorem 1 (The generating function) For any n > 1, we have 

A n (t) _ ^ k 

v 1 k>0 
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This theorem is a classic result in enumerative combinatorics, going back 
at least to work of MacMahon [H Chapter IV, §462]. See PU Section 4.5] 
for a modern treatment and generalizations. 

So, for example, the reader can check the n = 3 case gives the generating 
function for the sequence 1, 8, 27, . . ., of cubes: 

t + At + J =t + 8t 2 + 27t 3 + 64t 4 + 125t 5 + • • • . 
(1 — t) 

Another consequence of the balls in boxes approach is Worpitzky's iden- 
tity, which is explained in Knuth's book [131 Section 5.1.3]: 

l<i<n ^ ' ' 

This equation shows Eulerian numbers describe a kind of "change-of-basis" 
between binomial coefficients and nth powers. 
For example, when n = 4, k = 3, 

AiQ + A4 > 2 (Y) + ^' 3 u) +A ^a) = 1-15 + 11-5 + 11-1 + 1-0 = 81. 
To obtain Worpitzky's identity, we first recall from Equation ((IJ) that 



t l \^ (k + n — i 

(1 - t) n+l ~ 2-" V n 

Then from the left-hand side of Theorem (H we get: 

Mt) ( j \ 

n 

= ^2 An ^ ' 



i=0 

k + n — i 



i=0 k>0 

Tl 



n 



fc>0 \i=0 



But according to the right-hand side of Theorem [Q the coefficient of t k for 
this series is k n . Hence, (121) follows. 
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Next, we can derive a recurrence for Eulerian polynomials with Theorem 
[□and some calculus. Let F n (t) = J2k n t k = A n (t)/(1 — t) n+1 , and observe 
that 

F n {t) = Y,k n t k 

( nK-ijt) A' n ^(t) \ 
{l-t) n+1 

Comparing numerators yields the following well known identity for Eulerian 
polynomials: 

A n (t) = ntAn^it) + t{\ - ^A'^it). (3) 

From this equation, it is possible to deduce that the Eulerian polynomials 
factor completely over the reals, and moreover, the roots of A n (t) are all 
distinct, nonpositive numbers. According to Carlitz [1] this fact was first 
observed by Frobenius. 

By comparing the coefficient of f on the left and on the right of (El), we 
also get a handy recurrence for Eulerian numbers: 

A nA = iAn-tj + (n + 1 - z)yl n -i,i-i, (4) 

which resembles a weighted version of Pascal's recurrence for binomial coef- 
ficients. 

The reader should check this numeric recurrence against Table [1] to see 
how nicely it works. When drawing the numbers in a triangular array, we 
don't simply add the two numbers above to get the next entry, we add the 
appropriately weighted linear combination of the two. For example, the first 
few rows of such a triangle can be generated by hand as follows (so no need 
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to write out all 120 permutations in 65 to get the numbers A§ 




There is also a direct combinatorial proof of this recurrence. Imagine 
inserting the number n into a permutation of SVi-i- What happens if you 
insert it in a descent position? What happens if you insert it in an ascent 
position? See, for example, [121 Section 5.1.3]. 

There are many, many, more things that could be said at this point about 
Eulerian numbers. Some of these, including the row symmetry an astute 
reader may have noticed in Table [T], will be discussed at the end of the 
paper. 



Two-sided Eulerian numbers 

We are now going to extend the method of balls in boxes to study the "two- 
sided" Eulerian polynomial, 

n 

We think of A n (s,t) as the generating function for the joint distribution of 
left and right descents since A n ^j is the number of permutations w with i — 1 
left descents and j — 1 right descents, i.e., 

A n ,i,j = \{w E S n : des(w _1 ) = i — 1 and des(w) = j — 1}|. 

In Table [2] we see the two-sided Eulerian numbers for n < 8. 



A n (s,t) := sdeS[W 1)+1 t des ^ +1 



wes n 
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The two-sided Eulerian numbers were first studied by Carlitz, Roselle, 
and Scoville in 1966 [5], though rather than descents and inverse descents, 
they looked at the equivalent notion of "jumps" (ascents) and "readings" 
(inverse ascents). The results we prove here are all presented in Section 
7], and proved using a mixture of combinatorics and manipulatorics (i.e., 
manipulation of formulas using binomial identities and such). Here we take 
the balls in boxes approach to prove: 

• an analogue of Theorem [Q 

• a Worpitzky-like identity to compare with (jSJ), and 

• a recurrence relation that refines (J4J). 



Balls in boxes, revisited 



To prove these "two-sided" results, we will put balls in boxes another way. 
First, we introduce a two-dimensional analogue of barred permutations. 

Permutations in S n can be represented visually as an array of n indis- 
tinguishable balls so that no two balls lie in the same column or row. If 
w(i) = j, we put a ball in column i (from left to right) and row j (from 
bottom to top). For example, with w = 562143, we draw: 

6 
5 
4 
3 
2 
1 



j/i 



1 2 3 4 5 6 



The advantage of such an array is that w' 1 is easily seen in this picture. 
Rather than reading the heights of the balls from left to right to get w, we 
can read the column numbers of the balls from bottom to top to get w" 1 . 
So in the example above, w~ l = 436512. (Check this! Why is it generally 
true?) 

The analogue of a barred permutation, which we call a two-sided barred 
permutation, is any way of inserting both horizontal and vertical lines into 
the array of balls, with the requirement that there be at least one vertical 
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n = 6 
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Table 2: The two-sided Eulerian numbers [A n ,i,j]i<i,j< n , n= 1, . . . , 8 
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line between balls that form a descent in w, and at least one horizontal line 
between balls that form a descent in w' 1 . Other horizontal and vertical lines 
can be added arbitrarily. For example: 



(5) 



In a sense, this two-sided barred permutation corresponds to two ordinary 
barred permutations, one for w. ||56|2||14|||3, and one for w' 1 : |4|3||6|5|12|. 

Now fix an arrangement of balls corresponding to a permutation to. We 
will give a two-sided barred permutation the weight 

^(^of horizontal bars -f-lj^f^of vertical bars 

so the example seen in (jSJ) would contribute s s t 9 . Since the horizontal and 
vertical bars can be inserted independently of one another, we see that the 
generating function for the number of two-sided barred permutations corre- 
sponding to a fixed permutation w in S n is the product of the generating 
function for barred permutations for w' 1 (in the variable s) with the gen- 
erating function for barred permutations of w (in the variable t). Thus (pQ) 
gives: 

gdes(iu )+l£des(«;)+l 

(l-s) n+1 (l-t) n+1 " ^ 

Adding up (jSJ) over all permutations in S n we get the generating function for 
all two-sided barred permutations is 



A n (s,t) 



(1 - s^+^l-t)^ 1 ' 



(7) 



Now consider forming two-sided barred permutations with the bars first. 
Given k — 1 vertical bars and I — 1 horizontal bars, we get a k-by-l grid of 
boxes in which to place our balls, with the convention that if more than one 
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ball goes into a particular row or column, we arrange the balls diagonally 
from bottom left to top right. For example, the following arrangement of 7 
balls in a 5 x 4 grid of boxes yields the following two-sided barred permutation 
(with underlying permutation 1723465): 



© 














© 
















© 


© 


© 
© 




© 



< — > 



With n unlabeled balls and kl distinct boxes, this means there are a total 



of 



kl\\ _ fkl + n-V 
n J J \ n 

two-sided barred permutations with k — 1 vertical lines and I — 1 horizontal 
lines. (We are essentially choosing n of the kl boxes, with repetition allowed.) 
The generating function for the number of all two-sided barred permuations 
is thus 

'kl + 71- 1\ M 

n 



£ 

k,l>0 



S t 



and comparing with (j7|) yields the following. 

Theorem 2 (The two-sided generating function) For n > 1, we have 

W = >p (kl + n - 1\ k l 

Just as Worpitzky's identity can be derived from Theorem [T] and Equa- 
tion ([TJ, Theorem |5] and Equation ([T]) yields the following "Worpitzky-like" 
identity of binomial coefficients. 



kl+n-1 
n 



^1 A-n,ij 



k + n — i\ (I + n — j 
n J \ n 



(8) 
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To get dSJ), we use (Q]) and flU]) to see that for any permutation w E S n 
with des(u> -1 ) = z — 1, des(w) = j — 1, the generating function for its barred 
permutations is 

s l P sr-^ fk + n — i\ k fl + n — j\ l 



(1 - s) n+1 (l - t) n+1 ^\ n J ^ V n 
y / V / fc>0 v 7 z>o v 



n I \ n 

k,l>0 v 7 v ' 

There are Ai,i,j such permutations, so summing over w £ S n yields 
A n (s,t) . / sH* 



U = l v v 



;i - s) n+i (i - t) n+i f-± x n ' lJ V(i - s) n+1 (i - t) n+1 

i,j=l k,l>0 v 7 v 7 



E(EM* + rO('T')) A 



fc,Z>0 \tj=l 



as the generating function for all two-sided barred permutations for S n . But 
we already know there are ( fc ' + ™ _1 ) two-sided barred permutations of weight 
s k t l , and so we have (jSJ), as desired. 

We remark that this Worpitzky-like identity can also be proved (and con- 
siderably generalized) using the method of bipartite P-partitions, as follows 
from [121 Corollary 10]. 

Now from Theorem [2] we obtain a recurrence relation akin to ([3]), which 
makes for easy computation of the two-sided Eulerian polynomials. See [5j 
Equation 7.8]. We have: 

nA n {s,t) = (n 2 st + (n- 1)(1 -s)(l - t))A n ^{s,t) 

d d 
+ nst{\ - s)—A n ^(s,t) + nst(l-t)—A n ^(s,t) 

+ st (l_ s )(l_ t )_^_A ri _ 1 ( S ,t). (9) 

To obtain ©, we let F n (s,t) = A n {s,t)/{1 - s)" +1 (l - t) n+1 , and note 
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the following identity of binomial coefficients that we will need: 



n 



kl + n-1 
n 



n- (kl + n- 1)! 
n\(kl - 1)! 

(kl + n- 1) 



(JW + n-2)! 



(n-l)!(JW-l)! 



kl 



kl + n -2 
n — 1 



+ n 



kl + n -2 
n — 1 



Thus, we can see that 



nF n (s,t) = ^ n [ 

k,l>0 ^ 



k,l>0 
d 2 



kl + n-1 

n 

kl + n -2 
n — 1 



fc,i>0 



kl + n - 2 
n — 1 



st- 



F„_i(s,0 + («-l)^n-l(s,t). 



(10) 



Now, with a little calculus, we have 

nA n _!(s,t) 



;i - s ) n (i - t) n+1 1 (i - s) n (i - t) n 



+ 



n£A ra _i(g,t) 
(1 - s) n+1 (l - t) n+1 ' (1 - s) n (l - t) n 



n 2 yl n _i(s,t) 



+ 



+i 



nfA n _i(s,t) 



+ 



dsdt 



A n _i(s,t) 



(1 _ _ f)n (l-s)"(l-t)n 

After giving each term a common denominator, we compare numerators 
in (fTUj) to get as desired. 

By comparing coefficients on both sides of we also get the following 
four-term numeric recurrence: 

nAn,i,j = (ij + n- l)A n _i,ij 

+ [l-n + j(n + l- i)] + [1 - n + i(n + 1 - j)] Ai-mj-i 

+ [n - 1 + (n + 1 - i)(n + 1 - j))A n ^ ld ^. (11) 

The reader is invited to use this recurrence to obtain the first few arrays 
in Table EJ Is there a nice visual way to understand this recurrence? 
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Valley-hopping and 7-nonnegativity 

A sharp-eyed reader may have already noticed the well-known fact that the 
Eulerian numbers are symmetric for fixed n: 

-™-n,i ■^•n,7i+l— ii 

or in terms of Eulerian polynomials, A n {t) = t n+l A n (l/t). This symmetry 
can be explained by observing that if w has i descents it has n — l — i ascents. 

Similarly, one can find the following symmetries for two-sided Eulerian 
numbers in Table [2J 

A n ,i,j = ^n,j,ii ° r A n (s,t) = A n (t,s), (12) 
An,i,j = AvH-i-i.n+i-j, or A n (s, t) = (st) n+1 A n (l/s, l/t) , and (13) 
= A nin+ i_ J>+ i_i, or A n (s,t) = (st) n+1 A n (l/t,l/s). (14) 

It is a fun exercise to come up with combinatorial arguments to verify 
symmetries ffl~2l) and f|T3|) . (Hint: how does flipping a permutation array 
upside-down affect descents and inverse descents?) Symmetry ffT4l) follows 
from the first two. 

Another property exhibited by the Eulerian polynomials is unimodality, 
i.e., the Eulerian numbers in a given row increase up to a certain maximum 
and then decrease: 

Most good distributions satisfy this property of having the bulk of the mass 
is in the middle and the rest is spread out symmetrically. The canonical 
example is the binomial distribution, ( n ) with fixed n. Similarly, the two- 
sided Eulerian numbers appear to increase in the direction of the middle of 
the main diagonal, with a maximum at A n ^ \ n /2\ , \ n /2\ ■ However, things are 
somewhat delicate. For example, when i < j < n/2 it is not always true 
that A nyi+ ij > A n ^ t j < A n> jj_i, i.e., moving toward the main diagonal can 
sometimes lead to a smaller number. The first examples of this occur for 
n = 8, e.g., A 8 ^ 3 = 126 > 84 = A 8i2)2 , and A^ A = 1980 > 1773 = A 8>3j3 . 
See Tabled 

Both the symmetry and unimodality of the Eulerian polynomials follow 
from a stronger result, first proved in 1970 by Foata and Schiitzenberger [T0| 
Theoreme 5.6]. See also Carlitz and Scoville [6]. 
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Theorem 3 Forn > 1, there exist nonnegative integers 7^, i = 1, . . . , \n/2] , 
such that 

Mt)= E ln,t\l + t) n+1 - 2 \ 
l<i<\n/2] 

For example, when n = 4, 5 we have: 

A 4 (t) = t + lit 2 + lit 3 + t 4 

= t(l + t) 3 + 8t 2 (l+t) } 
A 5 (t) = t + 26t 2 + 66t 3 + 26t 4 + t 5 

= t(l + t) 4 + 22t 2 (l + t) 2 + 16t 3 . 

In other words, the polynomials A n (t) can be expressed as a positive sum 
of symmetric binomial terms with the same center of symmetry. Since the 
binomial distribution is symmetric and unimodal, so is the Eulerian distri- 
bution. One might say the Eulerian distribution is "super binomial" . Gessel 
has conjectured a similar expansion for two-sided Eulerian numbers that we 
will discuss shortly. 

First, we explain how Foata and Schutzenberger's result can be given 
a wonderfully visual proof via an action called valley- hopping. The argu- 
ment here has its roots in work of Foata and Strehl |TT] from 1974, was 
re-discovered by Shapiro, Woan, and Getu [15] in 1983, and was dusted off 
more recently (and applied in a wonderful variety of ways) by Branden [3] in 
2008. 

We define the valley- hopping action on permutations as follows. We vi- 
sualize permutations as arrays of balls again, but now we connect the dots to 
form a kind of mountain range. Some balls sit at peaks, others sit in valleys, 
and the rest are somewhere in between. If a ball is not at a peak or in a 
valley, it is free to jump straight across a valley to the nearest point on a 
slope at the same height. 

Valley hopping naturally partitions S n into equivalence classes according 
to whether one permutation can be obtained from another through a sequence 
of hops. For example, the permutation w = 863247159 would be drawn as 
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follows: 



t+i 




(15) 

There are 2 6 permutations in its equivalence class, formed by choosing which 
of the six free balls will be on the left sides of their respective valleys and 
which will be on the right. Write u ~ w if u can be obtained from w through 
valley-hopping. Notice that when a free ball is on the right side of a valley, 
it is not in a descent position, while if it is on the left side of a valley, it is 
in a descent position. (This is why we labeled the arcs in the picture with 
t + i.) Moreover, this property holds true regardless of the positions of the 
other free balls. 

As for the non-free balls, we know that peaks are always in descent posi- 
tions while valleys are never in descent positions. If a permutation has % — \ 
peaks, then it must have i valleys, and the remaining n + 1 — 2i balls are 
free. Thus, we can conclude that for a fixed w in S n with i — 1 peaks, 

'y ^des(w)+l _ _j_ 2i 

For example, the equivalence class for w = 863247159 would contribute 

t 2 (l + t) 6 . 

Since the union of all equivalence classes is S n , we see that the Eulerian 
polynomial is a sum of terms of the form t*(l + t) n+1_2j , proving Theorem 
|3j Moreover, the coefficient j n> i equals the number of distinct equivalence 
classes with i — 1 peaks. 

We now turn to the two-sided case. 
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It is possible to show that any polynomial in two variables (of degree n 
in each) satisfying symmetries (I12p and ( 113)) can be written uniquely in the 
basis 

{{st)\ s + ty(i + st) n+1 -i- 2i } < j+2t < n+1 . 

Thus, the two-sided Eulerian polynomials can be expressed in this basis, 
and Gessel has conjectured that such an expression is nonnegative. See [3j 
Conjecture 10.2]. 

Conjecture 1 (Gessel's conjecture) For n > 1, there exist nonnegative 
integers J n ,i,j , < i — j + 2i < n + 1, such that 

A n ( s ,t) = J2in,Asty(s + ty(i + st) n+1 -^- 2 \ 

For example, when n = 4, 5, we have 

A 4 (s, t) = st + 10(st) 2 + 10(sf) 3 + (st) 4 + s 2 t 3 + s 3 t 2 

= st(l + st) 3 + 7(st) 2 (l + st) + (st) 2 (s + t) 
A 5 (s, t) = st + 20(st) 2 + 54(st) 3 + 20(st) 4 + (st) 5 + 6sH 3 + 6s 3 t 2 + 6s 3 t 4 + 6s 4 t 3 

= st(l + st) 4 + 16(st) 2 (l + st) 2 + 16(st) 3 + 6(st) 2 (s + t)(l + st) 

In terms of the arrays [^4n,i,j] , we see: 

"10 

10 1 

1 10 

_ 1 

and 
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It is possible that one might use a "manipulatorics" approach to prove 
Conjecture [TJ (Perhaps an inductive proof using the recurrence in ©?) 
However, a more satisfying proof might be one that generalizes the valley- 
hopping proof of Theorem [3l 

If we look at both descents and inverse descents for the valley-hopping 
equivalence class of w = 12 ■■■n (i.e., the class with no peaks), we get a 
distribution of st(l + st) n ~ l . This is encouraging, but for the class of w = 
863247159 shown in (fT5l). we get 

s des(.-l) + l t des(n)+l = ^ + f) 2 ^ + ^ 

which is not even symmetric in s and t. So valley-hopping as presently done 
does not immediately give us a way to prove Gessel's conjecture. 

How should we partition S n so that we get groupings whose distribution 
of descents and inverse descents is given by (st) l (s + t) J '(l + st) n+1 ~^~ 2l l 

Generalization 

We finish by remarking that S n is an example of a finite reflection group, 
or Coxeter group. The notion of a descent can be generalized to any Coxeter 
group, and there is a "Coxeter-Eulerian" polynomial that enjoys many of the 
same properties of the classical Eulerian polynomial, including an analogue 
of Theorem [31 Moreover, this polynomial has topological meaning, as the 
"/i-polynomial" of something called the Coxeter complex. See [21 [71 [17] . 

The two-sided Eulerian polynomial generalizes to Coxeter groups as well, 
and seems to enjoy many of the same properties of A n (s,t). In particular, 
the analogue of Conjecture [H appears to hold in any finite Coxeter group. It 
would be interesting to have a general approach to the problem. 
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